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ABSTRACT 


Eqiuations  are  derived  for  detersinlng  the 
direction  cosines  of  a  plane  or  spherical  wave 
f^ont  hy  sissixadng  the  a2:*rival  tine  at  each 
■icroptione  in.  an  array  to  he  eun  Independent 
observation  and  reqyilrlng  that  the  sum  of  the 
squares  of  the  corrections  to  the  Irtdlvldu&l 
recordings  be  a  mini  mm. 

In  addlbloc,  the  effects  on  the  direction 
cosines  and/or  the  tine  arrival  errors  result¬ 
ing  from  considering  speed  of  sound  variations 
(profile)  vlhh  height  are  also  considered. 

In  particular,  five  examples  were  considered 
for  different  profiles,  none  of  which  were  extreme, 
and  resTilted  In  direction  angle  errors  as  large 
as  two  degrees  and  apparent  time  errors  on  the 
order  of  .04  second. 
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INTRODUCTION 


Various  publications  appear  in  the  literature  /l,  2,  ^  in  which 
:c]ethods  are  presented  for  detenaining  the  direction  (azimuth  and  elevation) 
of  propagation  of  a  sovmd  wave.  Only  a  few  attempts  are  made  to  correct 
for  timing  errors  end,  in  general,  no  attempt  is  made  to  correct  for  wind 
and  temperature  variations  within  the  array.  This  is  probably  becaiise  of 
the  assumption  that  the  variations  due  to  timing  inaccuracies  are  greater 
than  the  meteorological  variations.  Dean  /K7  recently  presented  a  least 
squares  method  for  determining  the  direction  cosines  of  the  wave  by  mini¬ 
mizing  the  effects  of  timing  errors. 

The  purpose  of  this  report  is  twofold:  (l)  to  present  a  general 
derivation  of  equations  that  can  be  used  for  determining  the .  dlridtihon 
cosines  of  a  plane  or  sphericeJ.  wave  front  when  considering  the  vertical 
meteorological  profile  over  the  array  and  mlsdLmizlng  the  time  errors;  and 
(2)  to  present  the  results  of  several  exanpleB  In  which  wind  and  tsmpcra- 
ture  profiles  over  the  array  are  considered. 


DISCUSSION 

A  least  squares  solution  and  a  linear  solution  are  derived  for 
determination  of  the  direction  cosines  of  a  plane  or  spherical  wave  travel¬ 
ing  within  a  sound  ranging  array  where  the  meteorologloal  parameters  vary 
as  a  function  of  height.  The  solution  is  based  on  but  not  restricted  to 
a  square  arrsy  of  four  microphones,  one  at  each  comer  of  the  array.  All 
computations  of  exaiaples  pmsented  are  referred  to  one  of  the  microphones 
of  the  array  as  a  base. 

Vertical  profiles  of  metf''— lata  (in  some  cases  hypothetical 
and  in  others  actual)  were  introduced  into  the  equations  and  the  time  of 
arrival  of  a  plane  wave  ITont  at  each  microphone  was  obtained.  Using 
these  times  and  restricted  meteorological  data  (surface  data  only),  an' 
Inverse  set  of  direction  angles  was  obtained.  The  difference  between 
the  two  sets  Is  considered  the  error  due  to  neglecting  the  total  vertical 
profiles .  (Examples  are  presented  which  reflect  the  order  of  magnitude 
of  these  errors,  but  unfortunately  no  data  were  available  which  would 
reflect  extreme  errors.) 

The  observation  equations  used  In  the  following  type  of  adjustment 
are  bsiscd  on  the  aeeumptlone  of  the  eo-oalled  "direction"  method  of 
geodesy  /^.  Each  time  recording  Is  treated  as  an  Independent  obser¬ 
vation,  and  the  sum  of  the  squares  of  the  corrections  to  the  individual 
recordings  is  to  be  made  a  minimum.  A  single  time,  however,  taken  by 
itself  determines  nothing,  for  If  each  of  the  recorded  times  bo  changed 
by  a  constsuit  amount,  the  new  set  has  the  same  significance  as  before. 

The  effect  Is  simply  a  change  In  zero  time,  which  is  purely  an  arbitrary 
matter.  If  a  set  of  time  corrections  for  an  array  has  been  deteznlnad 
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in  any  manner  and  the  mean  of  these  corrections  is  not  zero,  the  sum  of  the 
squares  of  uhese  corrections  can  always  be  diminished  by  subtracting  from 
each  correction  the  mean  of  all  the  corrections  so  that  the  algebraic  sum 
of  the  reduced  corrections  is  zero.  Hence  in  any  least  squares  adjust¬ 
ment  of  the  time  recordings  of  an  array,  the  algebraic  sura  of  the  time 
corrections  must  be  zero.  To  allow  for  this  change  of  time  base  orrfor 
a  constant  correction  to  all  recorded  times  in  an  array,  an  unknown  constant 
correction,  w,  enters  into  all  equations  expressing  the  results  of  obser¬ 
vations  at  an  array. 

As  is  usual  in  this  type  of  derivation  ,  the  observation  equation 
will  be  written 

+  dti  -t§+w-Vi=0  (1) 

where:  tj  is  the  actual  (observed)  time  of  recording  at  the  i^  station, 
t^  is  the  computed  time  of  arrival  at  the  i—  station, 
dt^  is  the  change  in  recording  time  at  the  i —  station  due  to 


a  email  change  in  the  direction  cosines  1,  m,  n  of  the  wave-front 
(salject  to  the  condition  that  1^  +  m^  ^  n^  =  l). 


w  is  a  constant  time  correction  for  all  stations  of  the  array. 


(uul  V|  i£.  the  resldvial  at  the  i —  station.  It  is  the  sum  of  the  squares 
oi  the  Vj^'s  which  is  to  be  made  a  minimum. 

:dnce  the  origin  of  time  is  arbitrary,  time  can  bo  measured  convene 
lently  from  any  recording,  and  for  this  station  the  arrival  time  and  the 
calculated  arrival  time  will  both  be  zero.  Also  for  this  station,  as 
can  be  seen  later,  v  =  w. 

The  basic  equation  connecting  time  and  direction  cosines  will  first 
te  developed  for  the  simplest  of  all  possible  cases,  that  of  a  plane 
wave  propagated  at  a  constant  velocity  in  a  still  atmosphere,  and  later 
modified  to  more  complicated  cases. 
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If  a  plane  wave  (unit  naraal  having  direction  cosines  1,  n;,  n) 
passes  Station  T  at  time  t^  -  t^  -  0  and  is  propagated  with  a  speed  Cj, 
Its  calculated  tlae  of  arrived  at  Station  I  will  he 

"  c  ^*1  ■'  ^f  "  ^f^“  “  ^f^®  ^f  ^^f  “ 

Here  X,  Y,  Z  are  station  coordinates.  Since  n  is  a  function  of  1  and  ■ 

(n  m  yfi  -  1^  -  ) ,  substituting  In  the  above  equation  leads  tc; 

ti  -  -i  (Xj^  X^)l  +  (Y^  -  Y.^)m  ±  (Z,  -  Z^)  (l  1®  -  nl®)*  .  (j) 


This  Is  not  linear  In  1  and  m.  Expanding  In  a  Taylor's  series  about  an 

assumed  set  of  dl^rectlon  cosines  (1  .  m  ^  n  ),  noting  that 

o  o  o 

n  -  n^  -  (l^/n^)dl  -  (M^/n^)da, 
and  (Just  to  reduce  the  nundser  of  terms)  setting 


(“if’o  ■  <*1  -  ^f>"c  *  (‘‘j  -  'f>“< 

■  '*1  -  -  hW 

(''if’o  ■  <*1  -  Xf’  ■  <==1  ■■  ^fV-c 


(4),(J4a)* 


gives 


‘i  ‘  J  ">if >0  »  Z  <'if >0"^  ^  Z  fir'e**- 


Hote  that  for  any  given  least- sq^iares  adjustment  ^®if^o 

vlIlL  be  numbers  depending  only  on  the  station  coordinates  and 
the  assumed  direction  cosines,  but  must  be  recomputed  for  each  Itera¬ 
tion.  Substituting  back  into  equation  (l)  gives 

Vj  -v.i  (r.,)^d.  .i  (I>,,)o  -  t;.  « 


*  For  the  set  of  "a"  equations  vtalch  begin  on  page  ^  ,  equation  (4a)  le 
the  sane  as  equation  (4). 


These  equations  contain  w  which  is  of  no  Importance  in  itself,  being 
simply  a  constant  correction  to  aJ.1  times.  Since  the  sum  of  the  v_,  's  is 
to  equal  zero  (n  =  number  of  stations  ) 


0  (7) 


i 


(8) 


Substituting  (8)  into  (6)  gives  the  reduced  observation  equation; 
v^  *•  L^dl  +  Mjdm  +  (9) 


vnere 


■“1  ""  c^^“if^o  “ 

^i  °  ^^If  ^o  '■  N  ^°lf  ^o^ 


(t^ 


(10) 


As  before,  K^,  are  sets  of  numbers,  one  triad  for  each  station, 

which  depend  only  on  the  initial  conditions  and  must  be  recomputed  for 
each  iteration.  From  here  on  the  usual  solution  scheme  is  used.  In 
Gaussian  notation  the  normal  equations  are 


[lL]dl  +  [LM]dm  =  -  [LK) 

,LM]di  +  [MMjdm  =  -  [MK)  (11) 


and  by  Cramer's  Pule: 

^  im  [NK]  -  [IX]  [MM] 

[LL]  [MM]  -  [IM]  [IM] 

^ .  m  . 

[LL]  [MM]  -  [IM]  [IM] 


(12) 


4 


I 


A  new  set  of  direction  cosine*  is  calc\ilated.  froa 


1 '  -  1  +  dl 
0  o 

B*  -  m  da 

o  o 

«  ±  [  1  -  (r)*  -  (■^)*]*.  (03) 


Equations  (4)  and.  (lO)  are  used  to  set  xip  nev  observation  equations 
(9)  and  the  process  repeated  until  dl  and  dm  are  negligible.  Then  v 
is  calculated  from  equation  (6)  auid  the  v^'s  from  equation  (6). 

To  derive  the  equations  to  take  into  account  a  wind  of  constant 
velocity,  7,  pcunUel  to  the  earth's  surface,  assuae  V  •*  S^l  +  S^J.  This 

vlU  reduce  the  difference  in  arrival  times  at  Stations  I  and  7  by 


^  (S^l  +  SyB)  since  t^  -  0  .  (la) 

Then  equation  (2)  becoaes 

[(X^  -  X^)l  +  -  Yj)b  +  (Z^  -  Zj.)n]  -  (~)(S^1  +  Sya).(a*) 

Solving  for  t^  and  substltxiting  for  n,  gives 

t. - i -  [(X  -X  )i  +  {r  -Y  )m  ±  (z  -z  )  ( 1  -  i»  -  a*)*}/3*) 

^c  +  Sl  +  Sm  ^  ^  ^  ^  ^  ^ 

X  y 


This  again  must  be  linearized.  Letting 


c  +  S  +  S  a 
X  o  y  o 


(14),  (iW)* 


*  For  the  set  "a"  equatleas  vhlch  begin  on  this  page,  equatlea  (l4a)  is 
the  saae  as  equation  (l4). 
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tken 


■ubstitatlng  (%)  In  cqaation  (1),  and  tllalnating  v,  va  again  gat 
aqioatlon  (9) 

-  L^dl  -f  4 


Trcm  hart  on  aqiaationa  (11),  (12),  and  (13)  earry  an  tka  aolutlon* 


Ska  MavBVtloa  of  a  ooaataat  vlad  yalaoity  la,  la  gaaaral,  aaraallatle 
alaea  it  la  a  fvaetlon  af  bath  ayaea  and  tlM.  lovayar,  vkatarar  tka  fane- 
tlona  and  8^  vara  af  tlaa  and  ayaea,  tkara  aaat  axiat  a  (V*)^^  ■ 

^^x^lf^  *  ooaataat  vlad  yalaoitjr  vkiek  would  hara  tka  aaoM 

6 
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effect  as  the  variable  wind  during  the  test.  Parentheses  and  subscripts 
have  been  added  to  indicate  that  these  are  tc  be  tb.e  beet  estimates  of 
the  "average"  wind  velocity  components  during  the  time  Ir'.ervrl  {t^.  t^) 


and  over  the  path  the  sound  ray  travelled  to  Station  T. 
considered  to  be  functions  of  hei^t  ("  )  c-nly,  then^  to  a 
tion 

s.-i  r” 

^  h  Jc 


If  £  and  S  are 
X  y 


S  (h)  dh 

X 


S*  = 

y 


j 


S  (h)  dh  . 

y 


(15) 


If  it  la  desired  to  modify  these  equati'«ns  with  correctlcns  for 
wa've-front  curv;  ture,  then  informatlcn  must  be  obtained  regarding  the 
change  in  directlGC  :;03lne5  ever  the  array.  The  simplest  ascumpticn 
would  he  that  tne  wave  front  spreads  out  from  a  point  In  space  (X.-.y'yZ). 
Then  the  direction  cosines  to  Station  I  wo;i.ld  be 


cosa^  =  1 

^  L 

cosp^  =  m^^ 
cos7^  = 


(X,-  X)/R^ 

(z,  -z)/r^ 


where 


-  [(X  -  X^)®  •*  (Y  -  +  (Z 

Setting 


"i  “  "f  ^^If 
®i  “  ”f  *  ^If 

"i  ■  “r  *  “^if 


then 


A1 


X^  “  X  -  X 


If  R. 


(16) 


(IT) 


(^8) 


(.19) 
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To  apply  this  Bodel  It  vlll  he  necessary  to  determine  either  an  R 
or  a  point  X,T,Z  and  calculate  the  ^  ^  ^^s  .  Then 

equation  (4)  vlll  be  evaluated  for  each  station  with  respect  to  the 
chosen  base  station.  This  vlll  result  In 

Olf'ol  -  <*1  -  -  ’'fX-o  ♦  * 

(^1  -  ♦  a.!,) 

■  <*1  -  ’'f’  -  (^1  -  ><1„  ♦  (41,) 

•''if’ol  ■  <’'l  -  *(>  *  <^1  -  ♦  Ao„) 

vhere  1,  m^,  n  are  nov  the  assumed  direction  cosines  at  the  base 
o  0  o 

station  (Station  F). 

Then  equation  (?a)  vlll  be  changed  to 


.  .  r^^lf^ol  '°if^ol/-.4  L, 

'i  ♦  “i  ■  L~s - 


<Vol 

'ol  ''01 

^1 


^  °01 


(6.)  ♦  [i^] 

ol 


vhere 


'ol  ■ '  *  <'*>if  <^o  ♦  ^if>  ♦  <';>if  <■.  ♦  '"ifX 


Equations  (10a)  become 

.  r^®lf^ol  ^^If^olf-  j  1  iVT^^lf^ol  ^°lfW 


'ol 


'ol 


(5h) 

(14b) 


(10b) 
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and  again  substituting  (5b)  £uad  (lOb)  Into  equation  (l)  we  get 

-  L^dl  +  M^dm  +  K^.  (9) 

After  solving  for  dl  and  dm  using  equations  (U)  and  (12)  va  now  have 

1,  -  1  +  Al,  -  +  dl 
1  o  If 


m .  “  m  +  An .  *  +  dm 
1  o  if 


For  1  *  f  this  gives 


If  .  lo  .  dl 


■  m  +  dm 
f  o 


f  \1  1  -  i;;  •• 


These  will  be  the  new  assumed  direction  cosines  to  be  put  back  Into 
equations  (4b). 

When  the  Iteration  has  reduced  dl  and  dm  to  a  negligible  value, 
1^,  m^,  and  n^  w511  be  the  direction  cosines  of  the  wave-front  at 

Station  F.  If  the  direction  cosines  at  another  station  are  desired, 
they  can  be  computed  using  equations  (l8). 


Other  assumptions  regarding  the  curvature  of  the  wave-front  can  be 
put  into  these  equations  by  redefining  ^If '  ^If  equations 

(l8)  and  (19)«  The  form  of  the  equations  following  will  not  be  changed. 


o 
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EXAMPLES 


In  the  presentation  that  follows  It  has  heen  asstaasd  that  the  wave 
ffont  when  at  the  first  microphone  is  plane,  but  changes  orientation  as 
it  passes  over  the  array.  This  results  because  the  speed  of  sound  over 
the  array  is  not  constant  with  respect  to  elevation.  Using  a  modifi¬ 
cation  of  equation  (2a),  t^  is  determined  for  three  microphones  relative 
to  the  first.  In  turn,  these  values  of  ti  along  with  a  "loccd."  speed  of 
sc jnd  are  used  in  equation  (2a) '  below  to  calculate  the  direction  cosines 
which  are  then  compeared  with  those  InltleLlly  specified. 

Equation  (2a)  is  rewritten  as 

(X^  -  X^)l  +  (Y^  -  Y^)m  +  (Z^  -  Z^)n  -  (C  +  S^l  +  S^m)!^  =  0  (2a)' 

where  C  is  the  speed  of  sovind  due  to  temperature  only  as  a  function  of 
height;  Sx  and  Sy  are  wind  components  in  the  X  and  Y  directions 
respectively,  as  a  function  of  height- 

Letting 

C(h)  »  c  +  c  h 
S„(h)  -  a  +  a  h 

X  12 

S  (h)  =  b  +  b  h, 

y  12* 

(2a) '  can  be  written  as 

IX.  +  mY,  +  nZ.  -  +  c  h  -i-  l(a  +  a  h)  +  in(b  +  b  h)  Idt  ••  0  (20) 

1  1  iJLi2  12  isj 

t 

o 

where 


t  -  0,  and  X,  -  Y,  -  Z,  -  0. 
o  f  f  f 

This  can  be  simplified  by  using  the  average  values  of  C(h),  S^(h), 


and  S^(h),  i.e. , 


c  h 

)dh  »  c  +  ^ 


C(h)  •=  h  f  (c  +  c  h 

—  1  r  ^ 

S  (h)  ■  r  /  h)dh  =  a  +  * 

X  iij^  1  s 


i  (h)  -  c  ^  +  b  h)dh  -  b 

y  «Jo  1  a  1 


1  2 

a  h 

+ 

1  2 
b  h 
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Then  equation  (20)  reduces  to 

e  fc 
IX 


eh.  ah.  b  h  ^  ' 

.  +  mY.  +  nZ  -  c  +  —  +  l(a  +  — )  +  m(b  +  ■  ;  t  = 


Combining  like  terns,  the  above  reduces  to 
a  h 


Xi  -  t^(a,  + 


i  1  r  ^  ^'1 1  °  ^1 

-)  4  n  y.  -  t  (b  4  — )  4nZ.  -  t  (c  4  -* — )  >  0. 

a  2"  ^  •“!  ^ 


(21) 


Equation  (21)  contains  fou”  luknovns,  1,  m,  n,  and  h  .  For  a  square 
array  of  microphones  as  shown  in  Figures  1  and  2,  h.^can  be  determined 
and  written  as 


d. 

i -  (22) 

Ctn(7  “  |)  4  Ctn(K  -  y) 

where 

=  L  Sin  B  (L  Is  the  length  of  the  side  of  the  square) 

d  -  L  Cos  B 
2 

d  =  l/2  Sin  (  X-  4  B)  =  L(Cos  B  4  Sin  B)  =  d  4  d 
3  ^  12 

B  =  tan  ^  {—)  (Figure  5). 

P 

If  ^ 

angle  B  =  angle  A  =  ■^,  then  d^  =  d  =  •^-  L  and  d^  =  L  '/2 . 

Substituting  the  value  of  h^  (where  1=1,2,  end  5)  into  equation 
(21)  reduces  it  to  three  unknown^  i,  m,  and  n.  Having  three  microphones 
to  work  with  will  give  three  equations  in  three  iinknowns,  i.e.. 


lA 

+ 

mB 

4 

nZ  - 

C 

=  0 

1 

1 

1 

1 

lA 

4 

mB 

4 

nZ  •• 

C 

=  0 

2 

2 

2 

S 

lA 

4 

mB 

4 

nZ  - 

c 

=  0 

3  3  3  3 


where 


=  ^1  “(a^  4^)  t^ 

Ci  =  4 


(25) 
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FIG.  2 
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tL.ao) 


FIO.  3 


A  Is  the  angle  the  ware  front  mkes  with  the  y-axls. 
B  Is  the  angle  the  vare  front  ■akes  vlth  the  x-axls. 


TAB  A 


or  A  -  -5-  -  B 


01  Is  norsal  to  the  vare  front 
and  Is  wilty. 


MAW  M  COSOI _ m  ^  A 


OM  is  a  projection  of  the  noiaal 
onto  tho  xy-plane. 


11^ 
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CONCLUSIONS 


The  equations  presented  in  this  report  are  readily  adaptable  for 
stvidles  of  both  plane  and  spherical  waves  for  different  meteorological 
profiles.  Tout  are  necessarily  more  complicated  than  most  conventional 
methods  /3/» 

Besults  of  five  selected  exaisples  (see  Appendix  A)  show  that  errors 
an  large  as  two  degrees  occur  in  the  direction  ein^es  of  a  plane  wave 
If  one  Ignores  the  total  vertical  meteorological  (wind  components  and 
temperature)  profile,  l.e.,  to  a  height  h,  as  shown  in  Figure  2,  within 
the  ejTTty.  Possibly  even  more  significant  is  the  finding  that  when 
assuming  a  plane  wave  and  neglecting  most  of  the  total  verticsLl  profile 
(only  the  lowest  level  meteorologlceLl  data  were  used),  there  were  .. 
apparent  time  arrival  errors  as  large  as  .04  second. 

Even  more  Important  is  the  fact  that  these  results  do  not  reflect 
the  extreme  errors,  since  extreme  profiles  were  not  considered. 
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APPENDIX  A 


Results  of  five  examples  are  presented  in  table  form. 

All  colvmms  will  be  denoted  by  the  number  of  the  column  in 
parenthesis  and  interpreted  as  follows: 

1.  Original  selected  direction  cosines  and  azimuth 
of  a  plane  wave. 

2.  See  Figure  1. 

3-  Heights  at  which  meteorological  data  are  available. 

k.  East-west  components  of  winds  for  heights  in  (3)j 
(East  is  plus). 

5-  North-south  components  of  wind  for  heights  in  (3); 
(North  Is  plus). 

6.  Speed  of  sound  as  a  function  of  temperature  only 
for  heights  in  (3). 

7.  Arrival  times  calculated  by  using  values  in 
columns  (1+),  (5),  and  (6)  in  equation  (21. a)- 

8.  Angles:  Using  of  (7)  and  Sy.o,  '^yo,  and  Cq 
of(^)^  (5)  and  (6)  in  equation  21. a  for  obtaining  angles. 

9.  Times:  Using  a.  p.  /  of  column  (l)  and  ''-q; 

SyQ,  Cq,  from  columns  (4),  (5),  (6)  in  equation  21.a. 

10.  Time  arrival  errors:  (9)  -  (7)* 

11.  Direction  cosine  errors:  (8)  -  (l). 
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